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introduction 


Semiconductor crystals such as Hfri-£4J* grown by unidirectional 
solidification Bridgmann method have shown compositional segregations 
in both the axial and radial directions (Lehoczky et . al..l980, 1981, 
1983). Due to the wide separation between the iiquidus and the solidus 
of its pseudobinary phase diagram (Lehoczky and Szofran 1981), there is 
a diffusion layer of higher HgTe content built up in the melt near the 
melt-solid interface which gives a solute concentration gradient in 
the axial direction. The value of effective diffusion coefficient 
calculated from fitting of the data to ID model varies with 
growth conditions (Szofran 1984). This indicates that the growth 
condition of the as not purely diffusion controlled. Because of 

the higher thermal conductivity in the melt than that in the crystal in 
the growth system, there is a thermal leakage through the fused silica 
crucible wall near the melt-solid interface. This gives a thermal 
gradient in the radial direction. Hart (1971). Thorpe, Hutt and Soulsbv 
(1969) have shown that under such condition a fluid will become 
convectiveiv unstable as a result of different dif fusivities of 
temperature ana solute, It is quite important to understand the effects 
of this thermosolute convection on the compositional segregation in 
both axial and radial directions in the unidirectionally solidified 
crystals under various gravitational directions. To reach this goal, we 
start with a simplified problem to study the effects of thermal-solutal 
convection on the temperature and solutal fields under various 
gravitational orientations. We begin by reviewing model governing 
equations . 


Governing Equations 


In this study we adopt the Boussinesq approximation; The equation 
of state takes the form that density is constant except that in the 
presence of the gravitational field a buoyancy force exists due to 
density variations which is caused by the temperature variation and 
concentration variation in the melt. 


Under the Boussinesq approximation and axial symmetric boundary 
conditions, the governing equations in cylindrical coordinates for 
incompressible fluid flow of the system are: 
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Scale the dimensional variable 


The equations can be nondimens ionalized by scaling the variables 
by a factor F, i e V = FV* Scaling length by ^C, velocity by V&C, 
time by pressure by nondimens ionalize temperature by 


Q= 


t-t b 


settma Tm and nondimens ionalize solute concentration by setting 

Q) . After the scaling and dropping all the *, the dimensionless 
ecuations become : 
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Where the thermal and sclutal Grashof numbers respectively , Grp, Grc 
are defined by: 


GrrjMmi 


Gr c =i^ 


o‘ and « z [10] 

The Prandtl number, Pr, and Schmidt number. Sc, are defined by 

Sc=tw«c and Pr=«^*T. [11] 

Compare the nondimensionalized equations with the FIDAP 
equations, we use the following inputs to the FIDAP for strongly 
couoled eauations 
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Quantity 

Set tine 

Values used 

Density 

1 

1 

Viscosity 

1 

i_ 

Specific Heat,^P 

Pr Pr= V /« 

0.233 

Conductivity 

Kj / K, or 1 

2 

Caoacitv, 

1 

1 

Diffusivity, D 

1/Sc S c= v /D 

0.0143 

Thermal Volume 



Expansion, Pr 

G r ? G r r=gPrAT R|/ \P 

Solutal Volume 



Expansion, ft 

*1/ V 

> 


These governing equations show that the flow characteristic are 
determined uniquely by Git, Grc, p r and Sc. These equations have been 
solved by the FIDAP program developed by Fluid Dynamics International, 
Inc. The boundary conditions on the velocity field are no slip at all 
wails. The boundary conditions on the solute field are constant at the 
top of the melt and satisfy segregation condition at the growth 
interface . 


Conclusions 


Preliminary simulation results for the input values listed above 
and Grc=0 with Gtt =0 reveal that CdTe compositional profile under ID 
diffusion controlled growth condition agree well with the result 
obtained by Han et . al. (Kan et. al. 1992). Fixed grid simulation for 
Grc=0 with Gtt = 10 4 has also been obtained. Results indicated that CdTe 
concentration profiles has been effected by convection due to 
horizontal thermal gradients . (Figs 2). Although a great effort has been 
applied, the steady state simulations for the effects of concentration 
profiles under deformed grids has never been converged. The planed 
studies will be continued by doing transient simulations. 



Fig . 1 . Gr c =0 , Gr T =0 . Fig. 2 Gr c =0 , Gr T= io< 
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